act_site = [1 2 4 7 8];
pas_site = [3 5 6 9];

ind = 4;

for site = 1:9
    for layer = 1:4
        if cell_cnt{site,layer} > 0
            for cell_id = 1:cell_cnt{site,layer}
                cur_fb_cell_p = [prps_resp_by_cell{site,layer,:,cell_id}];
                cur_pb_cell_p = [prpb_resp_by_cell{site,layer,:,cell_id}];
                mean_prps_resp_by_cell{site,layer,cell_id} = mean(cur_fb_cell_p)-1;
                mean_prpb_resp_by_cell{site,layer,cell_id} = mean(cur_pb_cell_p)-1;
            end
        end
    end
end

mean_prps_act(ind) = mean([mean_prps_resp_by_cell{act_site,:,:}]);
mean_prps_pas(ind) = mean([mean_prps_resp_by_cell{pas_site,:,:}]);

sem_prps_act(ind) = std([mean_prps_resp_by_cell{act_site,:,:}])/sqrt(996);
sem_prps_pas(ind) = std([mean_prps_resp_by_cell{pas_site,:,:}])/sqrt(764);

mean_prpb_act(ind) = mean([mean_prpb_resp_by_cell{act_site,:,:}]);
mean_prpb_pas(ind) = mean([mean_prpb_resp_by_cell{pas_site,:,:}]);

sem_prpb_act(ind) = std([mean_prpb_resp_by_cell{act_site,:,:}])/sqrt(996);
sem_prpb_pas(ind) = std([mean_prpb_resp_by_cell{pas_site,:,:}])/sqrt(764);